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Abstract 

Using the convex optimization method and Helstrom family of ensembles introduced 
in Ref. we have discussed optimal ambiguous discrimination in qubit systems. 
We have analyzed the problem of the optimal discrimination of N known quantum 
states and have obtained maximum success probability and optimal measurement for 
N known quantum states with equiprobable prior probabilities and equidistant from 
center of the Bloch ball, not all of which are on the one half of the Bloch ball and all 
of the conjugate states are pure. An exact solution has also been given for arbitrary 
three known quantum states. The given examples which use our method include: 1. 
Diagonal N mixed states; 2. N equiprobable states and equidistant from center of the 
Bloch ball which their corresponding Bloch vectors are inclined at the equal angle from 
z axis; 3. Three mirror-symmetric states; 4. States that have been prepared with equal 
prior probabilities on vertices of a Platonic solid. 

PACS Nos: 03.67.Hk, 03.65.Ta 
Keywords: minimum-error discrimination, success probability, measurement, 
POVM elements, Helstrom family of ensembles, convex optimization, conjugate 
states 



1 Introduction 

In quantum information, the problem of detecting information stored in the state of a quan- 
tum system is of fundamental interest. In the simplest case, two dimensional systems or 
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qubits can be used to store quantum information. We assume that a quantum system is 
prepared in a certain state that is drawn with known previous probability from a finite set of 
known possible states and we want to find the best possible measurement that can be used 
to determine the actual state of the quantum system. If the states are mutually orthogonal, 
then they can be distinguished perfectly. But because of the quantum interference, it is 
impossible to discriminate quantum states by measurement. There are two basic approaches 
to accomplish state discrimination. In one approach, which is called the minimum- error 
discrimination and is the ambiguous discrimination, measurement outcomes are not allowed 
to be inconclusive but there is the possibility that the state is identified incorrectly. In this 
case, the probability of successful discrimination is made maximum by the optimum mea- 
surement. In the other base approach, which is called optimum unambiguous discrimination, 
no error occurs, but there exists a measurement outcome which gives an inconclusive result 
where we fail to identify the state. In this approach, it is tried to minimize the failure prob- 
ability by appropriate measurements. The topic of quantum state discrimination was firmly 
established in the 1970s by the pioneering work of Helstrom [2J, who considered a minimum 
error discrimination of two known quantum states and unambiguous state discrimination 
was originally formulated and analyzed by Ivanovic, Dieks, and Peres [3j HJ E] in 1987. |6J In 
this paper, we deal with minimum- error discrimination discrimination. We will use convex 
optimization j7] as a tool for reaching this aim, which has many other different applications 
where some of them have been seen in previous papers P, El El [101 EEH CGI [131 E!- 

In the present work, we have used convex optimization for discrimination of known quan- 
tum states with the aid of Helstrom family of ensemble idea in qubit systems. Applied 
optimization problem is minimization of upper bound of optimal success probability. Min- 
imum upper bound is equal with maximum success probability because there exist POVM 
elements that are orthogonal to conjugate states (by definition, any conjugate state have 
same convex combination with its corresponding known state) PQ. In this case, the POVM 
elements give an optimal measurement to discriminate given states. It has been proved, at 
least two conjugate states can be pure. If conjugate state is pure, then its corresponding 
optimal measurement is orthogonal to it and if conjugate state is mixed, then its corre- 
sponding optimal measurement is zero operator. We have also shown, it is impossible for 
all of Lagrange multipliers associated with the inequality constraints to take on the value 
zero. Optimal measurement elements and optimal success probability have been obtained 
for discriminating : 1. N equiprobable quantum states, which distance between each of them 
and center of Bloch ball is equal and, not all of which lie on the one half of Bloch ball, and 
all of the conjugate states are pure; 2. Arbitrary three known quantum states; 3. Diagonal 
N mixed states; 4. N equiprobable quantum states which distance between each of them and 
center of Bloch ball is equal and their corresponding Bloch vectors are inclined at the equal 
angle from z axis. 

The organization of this paper is as follows. In Sec. II, in summary we illustrate the 
measurement operators of constituting POVM for minimum- error discrimination of quan- 
tum states and Helstrom family of ensembles as a strategy for carrying out of discrimination. 
Then in Sec. Ill with transformation of the problem format to optimization problem we ob- 
tain the KKT conditions [7] and problem formulation, and then by using them we determine 
optimal measurement operators and obtain optimal success probability for N equiprobable 
quantum states located at equal distance from center of Bloch ball, which all of the states 
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are not on the one half of the Bloch ball and all of the conjugate states are on the boundary 
of Bloch ball and then we solve the problem exactly for two arbitrary quantum states, three 
arbitrary quantum states and some examples in Sec. IV which comprise: 1. Diagonal N 
mixed states; 2. N quantum states with equal prior probability and equidistant from center 
of the Bloch ball which their corresponding Bloch vectors are inclined at the equal angle from 
z axis; 3. The symmetrical mirror states l^i), IV^) an d \ips) that subject to transformation 
|0) — > |0), |1) — ► — 11) change into \4>2), \ipi) and \ip3), respectively; 4. Quantum states that 
have been prepared with equal prior probabilities on vertices of a Platonic solid. Finally, a 
brief conclusion and two appendices have been provided. 



2 Ambiguous quantum state discrimination and Hel- 
strom family of ensembles 

We assume that a quantum system is prepared with some known prior probability, in some 
state chosen from a finite collection of given known conceivable states. We want to identify 
the actual state of the quantum system. A state discriminating measurement determines 
probabilistically what the actual state of the system belongs to the set of possible states. 
We also assume, the state space 5 is a convex set in a real vector space and an operator 
ej on s is defined by an affine functional ej(pi) = Tr(piMjMj) from s to [0, 1] that Mj is 

quantum measurement operator and p(j\i) = Tr(piMjMj) is the probability to infer from 
the measurement that the system is in the state pj if it has been prepared in a state pi. The 

state of the system after the measurement is MjplJ f J — . The measurement operators satisfy 

J Tr(piMj Mj) 1 J 

the completeness equation 



N 



Defining IT, = M^Mj, then Uj is a positive operator such that [15] 



= I (2) 

i=l 

and p(j\i) = Tr(piHj). Thus the set of operators IT,- are sufficient to determine the proba- 
bilities of the different measurement outcomes. The operators Ilj are known as the POVM 
elements associated with the measurement that are suitable for minimum- error discrimina- 
tion. The error probability is expressed as 

N N 

Perr = ^ ^ p i Tr(p i U j ) (3) 
t=l j=l,j^i 

Suppose we are given a state chosen from {pi}fLi with a prior probability distribution 
{Pi}iLi(pi > 0, J2iPi = !)• O ur § oa l i s to find an optimal measurement to maximize the suc- 
cess probability to discriminate the states. It is sufficient to consider an N- valued observable 
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{Ilj}^ from which we decide the state was in pi when obtaining the output i. The success 
probability is 

N 

P = l-P err = Y^p t Tr( Pi U l ). (4) 

t=i 

The maximal success probability P opt is caused by the best operators {IT}^. In order to use 
the minimum-error discrimination strategy, we have to determine the particular detection 
operators {n i }^ 1 that maximize the right-hand side of the equation (j3J under the constraint 
([2]). We shall use a useful family of ensembles which have been introduced in Ref. [I] and 
is later shown to be closely related to an optimal state discrimination strategy. A set of 
N-numbers {pi, pi] 1 — pi, 7i}|=i is called a weak Helstrom family (of ensembles) if there exist 
N-numbers of binary probability discriminations {pi, 1 — pi\f = i and states {Ti}f =l satisfying 
^ = p < 1 and 

Pi 

piPi + (l - Pi)n = pjPj + (l - Pj)Tj (5) 

for any i,j=l,...,N. 

We assume that a priori probability distribution satisfies Pi ^ 0, 1 in order to remove 
trivial cases, p and r, are called Helstrom ratio and conjugate state to pi, respectively. It 
has been proved that [TJ 

P opt < p. (6) 

An observable {IT}^ satisfies P opt = p if Tr{rjli) = for any i = 1, N. In this case, 
the observable {n i }^ 1 gives an optimal measurement to discriminates {pi\f =l and we call 
the family {p i: p{, 1 — pi, Ti\f =1 Helstrom family of ensembles [1]. 
With 

Pi= -(I + bi.a), Ti = ~(I + Ci .a) (7) 

expression (jSJ) can be written in terms of c, and bj which are the corresponding Bloch vectors 
to pi and Ti, respectively, as 

pibi + (1 - pi)Ci = pjbj + (1 - pj)cj. (8) 

3 Problem formulation 

3.1 The case of N quantum states 

In this paper, we have restricted ourselves to qubit systems. In future, our method will be 
used in qutrit systems. 

We shall find optimal success probabilities and optimal measurements for discrimination 
of states pi, i = 1, N which have been prepared with prior probabilities Pi. We will see 
that minimum Helstrom ratio equals optimal success probabilities. 

Our problem is 

to minimize p, 
subject to |cj| 2 — 1 < 0, i = l,...,N; 
£ibi + (1 - pi)ci - pibi - (1 - pi)ci = 0, i = l,...,N 
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that have been formulated as an optimization problem (see Appendix A). It follows that this 
problem has the Lagrangian 



N 



i=i 

N-l 



+ 2j u M-2(plblx + (1 - Pl)%l ~ Pi+lb(i+l) x ~ (1 - Pi+l)x i+1 ) 
i=l 

N-l 

+ ^ ^3i-l(Pl&ly + (1 - Pl)Vl ~ Pi+lb(i+l)y - (1 - 



i=l 
AT-1 

+ ^(pi&i* + (1 - - p i+ i%+i) 2 - (1 - (9) 

i=l 

where b; = (6^, and c* = (a;,, t/j, 

The partial derivative of the Lagrangian with respect to p and x^ y^, Zi, (1 < i < N) must 
vanish. Thus, the KKT conditions with respect to z/j = (^-5, ^-4, ^-3), (2 < i < N) are 

|ci| 2 -l<0, i = l,...,N; (10) 

Pibi + (1 -pi)ci -Pibi - (1 -p<)ci = 0, i = 2,...,iV; (11) 
Ai>0, i = l,...,JV; (12) 

l + ^^.( Cl - Ci ) = 0; (13) 

i=2 

N 

2Aid + (1-^)5^^ = 0; (14) 

1=2 

2A i c i -(l-p < )^ = 0, i = 2,...,iV; (15) 

Aidcil 2 - 1) = 0, i = l,...,N. (16) 

By the relations ffTB"]) . ffT4l and (IT51) we can conclude that it is impossible which A, = 0, 
i — 1, N. Also, The KKT conditions conclude there can be at least two of c,, 1 < i < N, 

SO |Cj| = 1. 

Now our aim will be to solve KKT conditions to find optimal POVM elements. From 
(USD, (HU) and (TTSD 



N \ r 



(17) 



and 

Ai I C,; 1 2 1 



P< 2 
t=i ^ 
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and by calculating c* from (jSj) and then substituting it into (ITTj) we arrive at 

Cj = , r j = l,...,N. (19) 

Zv£=l (p-pO 2 

When |cj| = 1, we choose its corresponding measurement operator orthogonal to t$. 
Thus, with the aid of (ITTj) and (fT8j) the POVM elements are found as 

^ = ^-\Xi)(Xjl j = l,...,N (20) 

where 

\Xj)(Xj\ = l(I-Cr&) (21) 

and when |cj| < 1, the state is mixed and ITj corresponding to Cj is considered zero operator 
(and therefore A, = 0) in order that TrfaHi) = is satisfied for alH = 1, N and condition 
of P opt = p is provided. The terms corresponding to all states pj's that conjugate states to 
them are mixed states do not have contributions to the sum in the relation 

If states have been prepared with equal prior probabilities, then ffTTl) . (IT81) and f[T9l 
become 

N 



^A iC , = 0; (22) 



i=i 

v 



£ A 'W 2 = %r ; (23) 



8=1 



c, = — -—z ; A7r^~~T' i = 1, AT, (24) 



D b 

(^-i)Ef=iAr^-i 

respectively, that we have defined D = Ei=i Ajbj. Therefore 



|c fc | 2 -| Cj f = , bZ " b ^ + 2D '^- b ^ , jlk = l,...,N. (25) 
J ' (A^P-I) 2 (iVp-l) 2 EiiAi 



As a special case we suppose, all of the states pi,...,pN are not on the one half of the Bloch 
ball (N > 4) and their corresponding Bloch vectors have equal length of b. We also suppose 
|c,;| = 1 for all of the vectors Cj, % = 1, N and then we can result D = 0. Therefore, 



\Xi){Xs\ = l{I+^- [ ), j = l,.-.,N (26) 

and 

P° pt = P=^ + b) (27) 

where we have used |c$| = 1 for some i. 

In the next two subsections, we precisely work out the maximum success probability and 
the optimal POVM elements for ambiguously discriminating between any two-states, with 
prior probabilities p±, P2 and among any three-states, with prior probabilities p±, P2, pz- 
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3.2 The case of two quantum states 

Although the case of two qubit states is studied [TH [17] we find instructive to see how the 
known solution follows from our method. 
Equations (THj) . (fl~5"l) and (TTHT) in this case are simply 

2A lCl + (1 - p x )v 2 = 0, (28) 

2A 2 c 2 -(l-p 2 )z7 2 = (29) 

and 

(1 - p x )(l - p 2 ) - 2Ai(l - pa) - 2A 2 (1 - px) = 0, (30) 

respectively. 

The considerations |ci| = 1, |c 2 | = 1 and equations ( l28l) . ( |29l) and ( |30l) can be used to 



drive 

A^l^i, A 2 ^ ,3!) 



and substituting these values into equation ffl9|) gives 



P2b 2 -pibx 

Cl = 2 p-l ' C2 = ~ Cl " ( ^ 

Using | Ci | = 1 and paying attention to p > pi and p > p 2 we obtain 

P op '=P =^(l + |p a b 2 -p 1 bi|). (33) 

Note that 

n 1 = l(/- P2b2 - pibl .q), n 2 = l(/ + P2b2 ' Plbl .a). (34) 

2 V 2p-l ; ' 2 V 2p-l ; V ; 

The minimum error probability p™ n is found to be 

P^ = \{l-\p 2 h 2 - Pl h 1 \). (35) 

It can be written as 

P ™ = ~(l-Tr\p 2 p 2 -pi Pl \) (36) 

which was originally found by Helstrom [16] where \Tr\p 2 p 2 — PiPx\ is the trace distance. 
It is obvious that, the minimum error probability is achieved when II x and II 2 are the 
projectors onto eigenstates of p 2 p 2 — p\p\ that belong to eigenvalues P2 ~ P1 — ^ 2b2 ^ Plbl and 
+ \ P 2b 2 - Plbl \ ^ respectively [xg_ 
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3.3 The case of three quantum states 

We now want to obtain the exact solution for discrimination of three arbitrary known mixed 
states. We place N=3 on the equations (TT4|) . (TT5T) and (|T8|) then we have 

2A 1 c 1 + (l-p 1 )(z7 2 + z7 3 ) = 0, (37) 

2A 2 c 2 - (1 -p 2 )V 2 = 0, (38) 
2A 3 c 3 - (l-p 3 )z7 3 = (39) 

and 

(l-ft)(l-&)(l-03)-2Ai(l-&)(l-^ 

(40) 

If now we make the assumption |ci| = 1, |c 2 | = 1, |c 3 | 7^ 1 then (139]) becomes z?3 = and 
fj37jl . (JHHD and gOJ are led back to (j2g]), (j29l and pOl). respectively. In analogy to the two 
states case, the results are given by 

A, = 1^1, A 2 = i^2, A 3 = 0; (41) 

_p 2 b 2 -pibi _pibi-p 3 b 3 (p -pi)(p 2 b 2 -pibi) 

— , C 2 — — Cjj C 3 — + —. — ; — , ) 

2p-p x -p 2 ' P-Pz {P-P3){2p-Pi -P2) 

P opt =P=\(Pi+P2+ b 2 b 2 - pibil); (43) 

1 p 2 b 2 -pibx 1 p 2 b 2 -pib! 

n i=o( J -^ •° r )> n 2 = -(/ + .a), n 3 = 0. (44) 

2 2p - pi - p 2 2 2p - pi - p 2 

The cases of |ci| 7^ 1, |c 2 | = 1, |c 3 | = 1 and |ci| = 1, |c 2 | 7^ 1, |c 3 | = 1 have similar results 
for optimal success probability and the optimal POVM elements. 

Now we consider |ci| = 1, |c 2 | = 1, |c 3 | = 1. Using [pjbj— Pjb,-] 2 = [(1— pj)cj — (1— pj)cj] 2 , 
we arrive at 

(p - Pl f + (p - Pj ) 2 - fabj - Pj h 3 ) 2 
Cl -° 3 ~ 2(p- Pi )(p- Pj ) ' ^- 1 ' 2 ' 3 " ^ 

Substituting f j4"5l) for Ci.c 2 , Ci.c 3 and c 2 .c 3 into (see Appendix B) 

( Cl .c 2 ) 2 + (c!.c 3 ) 2 + (c 2 .c 3 ) 2 = 2(c 1 .c 2 )(c 1 .c 3 )(c 2 .c 3 ) + 1, (46) 

we obtain 
that 

L = 4 ( -pip 2 + pip 3 + p 2 p 3 - pl ) I + 4 (pip 2 - pip 3 + p 2 p 3 - £» 2 ) J + 4 (pip 2 + pip 3 - p 2 p 3 ~p\)K 

+2IJ + 2IK + 2JK -I 2 -J 2 - K 2 , 
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M = 2[(p\p 2 + pipl - p\p 3 - p\p\ - p\p 3 - P2PI + 2pl) - piK - p 2 J + p 3 I]I 

+2[(-p\p 2 - pipl + pfp 3 + P\p\ - p\pz - P2V\ + 2pa) ~ Pi K + P2J - PzI]J 

+2[(-p\p 2 - p\p\ - p\p 3 - p\p\ + p\pz + p 2 pl + 2p\) + piK - p 2 J - p 3 I}K, 

N = \-p\p\ + pipl + p 2 2 pl -p\ + p\K + p\j - p\l\l 

+ \p\p\ ~ V\v\ + P\p\ -pt+pjK -p 2 2 J + pjl] J 

+ [p\pl + p\p\ - p\p\ -p\- p\K + p\j + p\l\K - UK, 

I = (pibi -p 2 b 2 ) 2 , J = (pibi -p 3 b 3 ) 2 , K = (p 2 b 2 - p 3 b 3 ) 2 . 

Furthermore, by substituting Ci.c 2 , Ci.c 3 and c 2 .c 3 from fH5|) into fIB-il) . (IB-iij) and (IB— . 
we can therefore give Ai, A2 and A3 in terms of the known quantities. Thus required optimal 
POVM elements are determined by means of ( |T9l) . ( 1201) and ( T2TT) . 

4 Examples 

1. We will consider the situation in which given states are 

Pi = ~(/ + h z a z ) = \[{l + 6i»)|0)(0| + (1 - 6«)|1)<1|], i = 1, -., N. (48) 

Then relation (1191) follows that the x- and y- components of vectors Cj, i = 1,...,N are 
zero, and thus there are only two vectors of Cj, i = 1, N which have length 1. While the 
last — 2 of POVM elements can be chosen zero operator, the relation f|T7|) leads to the 
requirement 

= (49) 

P - Pi P-P2 

Using the z-component of the condition (jSJ) and combining the relations f|49|) . ([TBI and (|20|) . 
we can obtain the following optimal results. 

n 1 = |o)(o|, n 2 = |i)(i| if Pl b lz < P2 b 2z (50) 

and 

n 1 = |i)(i|, n 2 = |o)(o| if Pl b lz > P2 b 2z (51) 

and 

P° pt = \(pi+P2 + \p2b z2 -pib z i\), (52) 

where b\ z and b 2z are components that maximize the right side of the relation (1521) over all 
b iz . 

Similar to the two state case, the operators II! and U 2 are the projectors onto eigenstates 
of p 2 p 2 -pipi. 



as 



2. Let us consider N mixed states pj, j = 1, N with the corresponding Bloch vectors 
bj = (frsin 8 cos tpj, 6sin 9 sin tpj, bcos 9), j = l,...,N (53) 
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with uniform prior probability distribution (pj = j?, j = 1, N). By (1221) we see 

N 

X>^ = 0, (54) 



which implies that vectors Ci,...,Cjv are in the z = plane. Furthermore relation (fill) 
concludes that bj — b,; is parallel with Cj — Cj and 

P=T7 + 4tTI "TT) *7^i- ( 55 ) 

iv iv I Cj — Cj | 

Therefore p is minimum if J | = 1, i — 1, N. 

Equation fl25|) concludes, D is orthogonal to all bj — bj, i,j = 1, N, i ^ j. Thus 

N 

D = (0,0,6cos#^A;). (56) 



i=l 



By applying ( 1241) . we can show 

b cos 6k — b, frsin^ 



iVp-1 = -jv-3l( cos ^' sin ^'°) (57) 
and note that | c< j = 1 and p > pi, i — 1, N, we can obtain 



p = -^(1 + 6sin6»). (5i 



The corresponding conjugate states to pj, j = 1, N are 



\h) = cos(^)|0> + sin(^)e^)|l>, j = 1, N. (59) 

By substituting (1551) in equation (TSUI) we obtain 

„ 4A,-(l + &sin6>), , , . . ^ 

n,-= J \ . fl — - X, Xi, J = l,...,iV. 60 

We must have 

| Xj ) = cos(^)|0) + sin(^)e^|l) (61) 

because Tr(Tj\xj)(Xj\) = m ust be satisfied for all j = 1, N. Hence P opt = p and |Xj)(x.7'|> 
j = 1, .., N are optimal measurement operator for discriminating the states pj, j = 1, N. 
These results accord with given results in example 4 of Ref. pp. 

3. For now we will consider the example in which the input states are given by 

= cos#|0) +sin0|l) (62) 

|^ 2 > =cos0|O) -sin0|l) (63) 
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I^s) = |0>, (64) 

with prior probabilities Pi,P2{= Pi) arid p 3 = 1 — 2p 1 . On the one hand, if |ci| = 1, |c 2 | = 1 
and |c3 1 1 then the relation fT4"3l) gives 

P°p* = Pl (l + sm26), (65) 

on the other hand if |ci| = 1, | C2 1 = 1 and | C3 1 = 1, we can use fj45|) to show 

(p — pi) 2 — p\{l — cos 49) 



Ci.c 2 ~- 

C1.C3 = c 2 .c 3 
and fjl6|) evidently leads to 



{P-Pi) 2 

P 2 ~ PiP — P + Pi(l — 2pi) cos26> 
(p- Pl )(p + 2 Pl - 1) 

2 1 + C!.C 2 
(C1.C3) = 



(66) 
(67) 

(68) 



and by substituting fl66|) and fIBTl) in fl68|) . the optimal success probability can be attained. 
It is 



pop 



(1 — 2pi) (pi sin 2 + 1 — 2pi — pi cos 2 < 



(69) 



1 — 2pi — pi cos 2 6* 

Therefore, there are two regimes depending on p± and which coincide at p[ = 2 + C os6>(si 1 n6>+cos e) - 
If pi > p[ (pi < p[), the optimal success probability is given by the relation fl69l) (the relation 

ABED) m- 



In Ref. |19j . the measurement strategy of minimum- error discrimination has been used 
for obtaining the relations ( 1651) and ( 1691 . 

4. In this example, we consider states pi, i = 1, N with equal prior probabilities, which 
form the vertices of a Platonic solid centered at the origin with length of edge a. By using 
relations (126]) and (1271) . we have given optimal success probability and the Bloch vectors of 
conjugate states. 
Pyramid: 



Cube: 



Octahedron: 



Dodecahedron: 



Icosahedron: 



^popt _ 1 ' 

gpopt _ 1 ' 

6P°p* - 1 ' 
20P°p* - 1' 



i= l,-,4; F* = i(l + J|a) 



2 = 1,.. .,8; P op * 
i = l,..., 6; P opt 



A H a 

8 V 2 ; 



V^«0 

6 V 2 ; 



/=1,...,20; P op * = — (1 + ^a) 

20 V 3 ' 



\2P°vt _ 1 



1,...,12; po P , = 1 (1 + VSW5 
12 l 2^ 



a) 



(70) 
(71) 
(72) 

(73) 
(74) 
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5 Conclusion 



Using the idea of Helstrom family, there is one method for ambiguous discrimination. In this 
method, Helstrom ratio is considered to be the cost function, subject to resulted constraints of 
Helstrom family of ensembles. If KKT conditions associated with the optimization problem 
are satisfied, minimum Helstrom ratio will be equal with maximum success probability. At 
least, two conjugate states to known states can be pure. Every optimal non-zero POVM 
element is orthogonal to its corresponding pure conjugate state and all of the optimal POVM 
elements corresponding to all mixed conjugate states are zero operators. It is not possible 
for all of Lagrange multipliers associated with the inequality constraints to take on the value 
zero. Our method has been restricted in qubit systems. 

Form of state space is not a sphere for qutrit systems and the located states in boundary 
of state space are not necessarily pure, therefore using this way for qutrit systems seems 
different and it will be investigated in the future. 



Appendix A 

A summary of convex optimization 

An optimization problem [7] has the standard form: 

minimize /q (x) , 

subject to fi(x) ^ 0, i = l,...,m, 

hi(x) =0, i = l, ...,p, (A-i) 

where the vector x = (xi,...,x n ) is called the optimization variable and the function fo : 
R™ — > R the cost function. The inequalities fi(x) < are called inequality constraints, and 
the equations hi(x) = are called the equality constraints. 
The Lagrangian L : R™ x R m x R p -> R is 

m p 

L(X, X, V) = f (x) + 22 ^ifi( X ) + U i h i( X )- ( A_ii ) 

i=l i=l 

We refer to \ as the Lagrange multiplier associated with the ith inequality constraint 
fi(x) > 0; similarly we refer to vi as the Lagrange multiplier associated with the ith equality 
constraint h^x). The vectors A and v are (Ai, A m ) and {i>\, u p ), respectively. 

The dual function g : R m+P — ► R is defined as the minimum value of the Lagrangian 
over x that is written as 

g(X, v) = inf xeD L(x, A, u) (A-iii) 

The dual function yields lower bounds on the optimal value p* of the problem (lA-ij) : for any 
A y and any v we have 

g(X,u)^p\ (A-iv) 
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A natural question is: what is the best lower bound that can be obtained from the La- 
grange dual function? This leads to the optimization problem 

maximize g(X, v) 

subject to A h (A-v) 

This problem is called the Lagrange dual problem associated with the problem (lA-il) . The 
problem flA-il) is sometimes called the primal problem. We refer to (A, v) as dual optimal if 
they are optimal for the problem flA- vj) . 

The optimal value of the Lagrange dual problem, which we denote d*, is, by definition, 
the best lower bound on p* that can be obtained from the dual function. We have 

d* ^ p* (A-vi) 

We refer to the difference p* — d* as the optimal duality gap of the original problem (IA-il) . 

We now assume that functions f , f m , hi, h p are differentiable. Let x* and (A*, u*) 
be any primal and dual optimal points with zero duality gap. Then we have 

fi(x*)<Q, i = l,...,m 

hi(x*) = 0, i = l, ...,p 
A* > 0, % — 1, m 

m p 

V/o(ar*) + E A * V ^) + E WW*?) = 

i=l i=l 

Xtfi(x*) = 0, i = l,...,m (A-vii) 

which are called the Karush-Kuhn- Tucker (KKT) conditions. The condition X*fi(x*) = 
0, % = l,...,m is known as complementary slackness; it holds for any primal optimal 
and any dual optimal (A*, u*) (when duality gap is zero). 

The converse holds, if the primal problem is convex. In other words, if x, A, v are any 
points that satisfy the KKT conditions and fa are convex and hi are affine, then x and (A, u) 
are primal and dual optimal, with zero duality gap. 



Appendix B 

Proof of (1461) 

By taking dot products of (13"7|) by c x and c 2 , and then using them and equations 
and ([TBI , the Lagrange multipliers associated with inequality constraints can be expressed 
as 

x = (1 -Pi)[(ci.c 2 )(c 2 .c 3 ) - (C1.C3)] 

1 2[1 + ( Cl .c 2 ) - (d-ca) - (c 2 .c 3 )][l - ( Cl .c 2 )]' 1 " ) 
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(1 -p 2 )[(c 1 .c 2 )(c 1 .c 3 ) - (c 2 .c 3 )] 

(B-nJ 



2[1 + (ci.c 2 ) - (ci.c 3 ) - (c 2 .c 3 )][l - (ci.c 2 ) 



As = r (1^1I1±%^ (B-iii) 
2[l + (c 1 .c 2 )-(c 1 .c 3 )-(c 2 .c 3 )] 1 ; 



Taking dot products of fl37j) by c 2 and c 3 and in the same manner, we find 

(i-pOK^.c,) 2 -!] 



Ai 



A 



2[2(ci.c 2 )(ci.c 3 )(c 2 .c 3 ) - (ci.c 3 )(c 2 .c 3 ) - (ci.c 2 )(c 2 .c 3 ) - (ci.c 2 ) 2 - (ci.c 3 ) 2 + (ci.c 2 ) + (C1.C3)] 

(B-iv) 

(1 -p 2 )[(ci.c 2 ) - (ci.c 3 )(c 2 .c 3 )] 



2[2(C 1 .C 2 )(C 1 .C 3 )(C 2 .C 3 ) - (Ci.C 2 )(c 2 .C 3 ) - (Ci.C 3 )(c 2 .C 3 ) - (C!.C 2 ) 2 - (Cx.C 3 ) 2 + (Ci.C 2 ) + (C1.C3)] 

(B-v) 

A = (1 -p 3 )[(c!.C 3 ) - (Ci.C 2 )(c 2 .C 3 )] 

3 2[2(Ci.C 2 )(Ci.C 3 )(c 2 .C 3 ) - (Ci.C 2 )(c 2 .C 3 ) - (Ci.C 3 )(c 2 .C 3 ) - (Ci.C 2 ) 2 - (Ci.C 3 ) 2 + (ci.C 2 ) + (C1.C3)] 

(B-vi) 

These equations lead to the equation ( 1461) . 
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